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Abstract 

We investigate the role of thermal instability, arising from radiative cooling of an 
optically thin, dusty plasma, by linear stability analysis. The corresponding isobaric 
stability condition for condensation mode is found to be modified significantly by the 
degree of ionization and concentration of finite sized, relatively heavy, and negatively 
charged dust particles. It has been shown that though the fundamental wave mode 
is similar in nature to the one in absence of dust particles, a new dust-wave mode 
can propagate in such a plasma. It is conjectured that the presence of negatively 
charged dust particles may considerably affect the stability of various astrophysical 
structures against thermal instability, which can not be explained with the help of 
gravitational instability. 

Key words: thermal instability, dusty plasma 
PACS: 53.20.Ex, 52.30.-g, 98.58.-w 



The formation and existence of a large number of astrophysical structures 
such as interstellar clouds, solar prominence, localized structures in planetary 
nebulae etc. can not be explained by means of conventional gravitational in- 
stability [1]. However, it has been argued that thermal instability may be a 
reasonably good candidate, with can accelerate condensation, giving rise to 
localized structures which grow in density by loosing heat, mainly through 
radiation [1,2]. The first comprehensive analysis of thermal instability in a 
diffuse interstellar gas is first given by Field [2]. Subsequently, many authors 
considered the process of thermal instability in different circumstances e.g. 
in hydrogen plasmas [3], to explain formation of clumpy gas clouds in two 
phase medium [4], and thermal instabilities in photoioized interstellar gas [5]. 
Recently, Birk [7] has considered thermal condensations in a weakly ionized 
hydrogen plasma. 
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In recent years, however, it has been reahzed that heavier dust particles consti- 
tutes an ubiquitous and important component of many astrophysical plasmas. 
[7-9] including interstellar clouds, stellar and planetary atmosphere, planetary 
nebulae etc. It is, thus of particular interest to see the effect of dust particles 
on the role of thermal instability in astrophysical situations. One of the obvi- 
ous facts of introducing dust as a component of plasma is the propagation of 
very low frequency compressional modes [10]. We, in this paper, investigate 
thermal instabilities in a multicomponent plasma with electrons, ions, and 
dust particles, assuming that a ^ d ^ Xd, where a is the average dust size, d 
is the average inter-dust distance, and is the Debye length. We consider the 
dust particles to be negatively charged and cold, whereas the ions can loose 
energy through radiation by de-excitation. We further assume that electrons 
are isothermal. We introduce a heat-loss function [2] exclusive of thermal con- 
duction, which takes care of the radiative loss. For simplicity, we do not take 
into account the dust-charge fluctuation, which usually has a stabilizing ef- 
fect [6]. The resultant equations are two-fluid magnetohydrodynamics (MHD) 
equations with Boltzmann distribution for the electrons. The dust acoustic 
(DA) and dust-ion acoustic (DIA) modes [10] are ruled out by virtue of 'cold 
dust' assumption. We show that the isobaric condition for thermal condensa- 
tion mode is significantly modified by the inclusion of dust dynamics. We also 
find the existence of a new wave mode, solely because of the presence of the 
dust particles. 



We begin by writing out the equations for a unbounded, weakly collisional 
dusty plasma with a heat-loss function T), which represents the rate of 
net loss of energy per unit mass, through radiation. We consider the dust 
particles to be negatively charged and cold and the electrons obey Boltzmann 
relation. We assume here that that ion-neutral collisional cross-section is of the 
same order as that of the electron-neutral collisions, so that the electrons can 
be treated as Boltzmann distribution if {Ti/Te){vTe/vTi) ^ 1, and in practice 
this condition is satisfied in the cases we are going to consider in this work. 
The equations are. 
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where Xi is the ion thermal conductivity, n are ionization and recombination 
frequencies, and the other symbols have their usual meanings. The subscripts 
e, i, d respectively refer to the electron, ion, and dust populations. Equa- 
tions (1) through (6) are the continuity and momentum equations for dusts, 
Boltzmann relation for the electrons, and continuity, momentum, and energy 
equations for the ions. The set of equations is closed by the quasineutrality 
condition, 

ZiUi = rie + ZdUd, (7) 

where represents the ion and dust charge numbers, respectively. Note 
that we have taken into account the possibility of multiply ionized ions. The 

ion pressure is given by pi = riiTi (Tj are expressed in energy units) and we 
have assumed that F, the ratio of specific heats for the ions to be 5/3. The 
recombination loss term in the ion-continuity equation represents the electron- 
ion recombination on the surface of the dust particles [11] and vr oc Trr'^ndUicl, 
where is the dust radius and Cg — {Te/rriiY^'^ is the sound speed. 

We consider now a small perturbation of the form exp(— ia;t + ik ■ r). The 
linearized expression for electron, dust, and ion densities are given by 



where we have used the linearized form of the Boltzmann relation for the elec- 
trons Eq.(3), dust-continuity equation Eq.(l), and the quasineutrality condi- 
tion Eq.(7). Wc have used ujd = kcd and = ZdTeo/ndo- The quantities with 
subscript '0' and '1' are equilibrium and perturbed quantities, respectively 
and '~' represents normalized perturbed quantities. We have normalized the 
perturbed densities of the species by their respective equilibrium values and 
perturbed electrostatic potential energy e0i by equilibrium electron thermal 
energy, = e^i/Teo. The quantities Sj = fijo/riio are the ratios of equilibrium 
electron and dust densities to that of the ions for j — e,d. 

The linearized ion-continuity and momentum equations are given by 

- iuui + i{k ■ Vii) = -uii{nd + Ui) + lyiUeSe, (11) 
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where cui — kci is the ion-sound frequency, q = (Tjo/mj)^/^ is the ion-thermal 
speed, r = Z.iTpX)/TiQ, and pi = Pn/pio is the normahzed perturbed ion pres- 
sure. FinaUy we write down the hnearized energy equation for the ions as 
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where we have defined Tj = Tn/Tgo, ujp = kpCi, and ujt — {kx + k^/kK)ci. 
Following Field [2], we define the wave numbers kp^T,K as 
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The first two are the wave numbers of isothermal and isochoric perturbations, 
respectively and the third one is the reciprocal of mean free path of the heat- 
conducting particles. It has been assumed throughout the derivation that the 
perturbed ion thermal conductivity Xi — ^ and the equilibrium heat-loss func- 
tion jC{po — miriio^Tio) — 0, where Cp^r represents {dC/dT)p and {dC/dp)T, 
evaluated for the equilibrium state. 

Using Eqs.(8)-(13), the dispersion relation for this radiating dusty plasma can 
now be written as 
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where we have used the relation pi — Ui^Ti. 

Before numerically solving the dispersion relation Eq.(15), however, it is in- 
structive to analyze Eq.(15) under certain limits. In the limit of coUisionless 
electron-ion plasma with no dust particles and thermal radiation and zero 
heat-conductivity, Eq.(15) reduces to 
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which is just the ion-sound wave modified by the ion-charge number. We now 
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look at the expansion of Eq.(15) [2], in the hmit of weak coUisions {i'lr ^ 
and vanishing thermal perturbation i.e. oUp^T 0, which, for the condensation 
mode, is given by 
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where we have substituted uo = in and neglected uod- The characteristic velocity 
of the modes, we are interested, is q and in a plasma with cold dusts {T^ -C 
Tj^e), as this one, Ci/cd ~> 1. For example, in interstellar clouds, rii ~ 10~^, 
Ud ~ 10-^ T^^Tir^ 10 °K and with Zd ~ 10^ and ~ lO^m^, q/q ~ 10^ 
so that wc can safely assume cUd <^ To the leading order, we have the 
radiation condensation mode if. 



UJp + ZiUj > Ur 
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which reduces, in the coUisionless limit to 

Z 
1 



cup> {1 + ^t]ujt. (19) 



Note that for zero thermal-conductivity, ojt — kxCi and the above conditions 
are independent of k. The conditions given by inequalities (18) and (19) are 
similar to the isobaric instability condition given by Field [2], modified by 
collisions and dust particles. We observe that the presence of a significant 
dust population (smaller 5e) considerably modifies the instability criteria for 
the condensation mode and has a stabilizing effect. It is interesting to note 
that the degree of ionization has also a similar effect on the stability of the 
condensation mode which becomes important in certain situations. For exam- 
ple, in condensations observed in planetary nebulae, the thermal instability 
is strongly associated with collisions of electrons with N+ and N++ ions [2]. 
It should, however, be noted that the conditions (18) and (19) do not reduce 
to one given by Field even in absence of dust particles, which is because our 
consideration of two-fluid approximation. 

The inequality (18) also implies that the growth rate of the condensation 
mode is independent of the wave number in the short wavelength limit, when 
the thermal-conductivity is zero in opposite to a thermally conducting plasma 
(finite kx), where there is no condensation for sufficiently larger k. Plots of 
the normalized growth rate (a) = cu/ujp) versus normalized wave number {k — 
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Fig. 1. Normalized growth rate of the colHsionless condensation mode vs. wave num- 
ber. The sohd curves are with dust particles whereas the dashed ones are without 
any dust. The topmost curves are for k]^^ = kp/kx = and the others for fe^^ = 0.01 
and 0.1, respectively. 

k/kp) for the condensation mode are shown in Fig.l, in absence of coUisions. 
The growth rate of the mode is definitely smaller in presence of dust particles, 
as wc have guessed, but note the lengthening of the cut-off wave number. 
In fact, the knees of the lower curves with finite k^ prove to a bifurcation 
point in a coUisional plasma, as we shall see shortly. The parameters for Fig.l 
are kr = kr/kp = 0.2, = 1, = 100, T,o = T^o = 10 °K, md/rrii = 10^, 
rijo = 10~^ cm~^, and nao — 10~^ cm~^. These values are typical for interstellar 
clouds [12]. 

An approximation similar to Eq.(17) for coUisionless wave modes is given by 




The instability criteria for a growing wave mode is thus {uip + 2/2>ujt) < 0, 
which is just the iscntropic instability criteria, given by Field [2]. We can 
thus argue that the principal wave mode is unaffected by the presence of dust 
particles, but for a decrease in magnitude, dictated by the term (S+SZjr/^e)"^, 
which proportional to 5e/Zf. However, there is one more wave mode which 
is not apparent from the above analysis. We can detect the latter only by 
appropriately solving Eq.(15), numerically. In Fig.2, we show the the growing- 
wave modes in a collisionlcss plasma with considerable dust presence as in case 
of a planetary nebula [13]. There is a purely dust-wave mode (second figure 
in Fig.2) which is otherwise absent in a pure electron-ion plasma. The growth 
rate of the mode is, however, much smaller. Notice the cut-offs in all these 
curves. These figures are with kr — —2, k]^ — 0.01, Zi — 2, Z^ — 100 and 
rido = 10~^ cm~^, all other parameters being same as in Fig.l. The value of kx 
is typical in planetary nebulae. We have chosen a large negative value for kx 
so as to drive the wave mode unstable (see Eq.(20)). Note that a condensation 
mode will always accompany these wave modes (not shown) as the condition 
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Fig. 2. Normalized growth rates vs. wave number for the condensation modes in 
absence of coUisions. The first figure shows the usual wave mode present irrespective 
of the dust particles. However, the growth rate in presence of dust particles is 
considerably smaller (solid curve). The second figure shows the dust- wave mode 
which is absent in a pure electron-ion plasma. 
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Fig. 3. Growth rates in presence of collisions with kx = 0.2, k]^ = 0.01, un/up = 0.2, 
ui/ojp = 0, and = 10~^cm~^. Other parameters are same as in Figs.l and 2. 

for condensation instability is readily satisfied in this case (kr < 0). 

In Fig. 3, we show the growth rates of modes with collisional effects. The solid 
curves represent condensation modes and the dashed one represents the wave 
mode. We see a merger of the condensation modes to a growing wave mode 
at a bifurcation point near k ~ 4.0. 

We should emphasize that we do not have an explicit expression for the in- 
terstellar heat-loss function i2(p, T) in general, because of an wide variety 
of compositions of material in different interstellar regions. Besides, several 
heating and cooling mechanisms may be operative simultaneously and their 
dependence on local density and temperature may not be known exactly. Del- 
garno and McCray [14], in their classic review, have mentioned about differ- 
ent cooling processes in HI regions. A relatively concise expression for the 
heat-loss function in photo-ionized interstellar gases with high metallicity is 
given by CorbeUi and Ferrara [5]. In the simplest case, for radiative cool- 
ing, the heat-loss function can be assumed to be proportional to p"'p"^, where 
n = 5/2, m = —1/2 for free-bound transitions and n = 3/2, m = +1/2 for 
free-free transitions [3,14,15]. However, the exact form the heat-loss function 
is not important at this point. As a prototype example, we assume the heat- 
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Fig. 4. Range of instability in terms of temperature with Sg- The value of 
T = 4.0, 4.5, and 5.0, respectively for the upper, middle, and lower curves. 

loss function used by Field [2] for thermal condensations in planetary nebulae, 
which is represented by 

where T is the mean excess ionization temperature of N+ ions. The tempera- 
tures are given in the units of Tn+ = 21800 °K, the mean excitation tempera- 
ture of the N"*" level. Applying the instability condition (19) in the coUisionless 
case, we have, thermal condensation if, 

r_(f,c)<r<r+(f,c), (22) 

where ( — {1 + 1/Se) and 

T±{T, C) = [2C + A;kT(2 + () ± {[2C + ^xT(2 + (23) 
-8^^fC(2 + 3C)}'/'](4 + 6C)-^ 

For simplicity, we have assumed T^q = Tiq and Zi = 1. The above expression 
of T± reduces to usual one when C = 1 (single- fluid approximation). For kx — 
3/2, Zd = 100, Hio = 10-3 cm-3, n^o = lO'^cm'^, and T = 4.5, we have T± = 
1.3293 - 2.542, which in absence of dust becomes 0.7248 - 3.7254 [2]. So, the 
presence of dust provides a stabilizing effect against radiative condensations. 
However, the range of instability shrinks to zero at T = 3.902, much higher 
than 1.974, which is for the dust-less case. This implies that for a planetary 
nebula with a considerable amount of dust particles, to be stable against the 
radiative condensation mode, it has to be at a much higher temperature. In 
Fig.4, we show the range of unstable temperature with Sg. Note that a lower 
value of Se corresponds to higher dust concentration. 

To conclude, we have found that the instability criterion for the thermal con- 
densation mode is significantly affected by the presence of dust particles, which 
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can severely alter the stability regimes of many astrophysical condensations 
where a significant number of dust particles are present. The fundamental 
wave mode is similar in both cases with and without dust particles, though 
we now have another wave mode, which is due to the presence of dust parti- 
cles. Finally, we have shown that the conventional stability regimes of many 
astrophysical condensations may be considerably changed. 
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